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ABSTRACT 

Weak gravitational lensing provides a sensitive probe of cosmology by measuring the 
mass distribution and the geometry of the low redshift universe. We show how an 
all-sky weak lensing tomographic survey can jointly constrain different sets of cosmo- 
logical parameters describing dark energy, massive neutrinos (hot dark matter), and 
the primordial power spectrum. In order to put all sectors on an equal footing, we in- 
troduce a new parameter /3, the second order running spectral index. Using the Fisher 
matrix formalism with and without CMB priors, we examine how the constraints vary 
as the parameter set is enlarged. We find that weak lensing with CMB priors provides 
robust constraints on dark energy parameters and can simultaneously provide strong 
constraints on all three sectors. We find that the dark energy sector is largely insen- 
sitive to the inclusion of the other cosmological sectors. Implications for the planning 
of future surveys are discussed. 
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1 INTRODUCTION 

In the last few decades, a wealth of c osmological data 
(from large scale structur e ([P eacock 2005, 2dF) ; the cosmic 
microwave b ackground (iKomatsu et all |2009|. WMAP51; 



supernovae 
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dAstier et all 120061. S NLS; iMiknaitis et all 
Kowalski et alj 12008t h weak lensing 



I Schrabback et al ] |2009l )1 las revolutionised our vision of 
the Universe. In this concordance cosmology, initial quan- 
tum fluctuations are believed to have seeded dark matter 
perturbations in which the Large Scale Structure we ob- 
serve today has formed. Within this concordance model the 
Universe is composed only of a small proportion of baryons 
(4%), the rest being dark matter (25%, which can be hot or 
cold) and dark energy. 

One of the main challenges today is to understand the 
nature of the mysterious dark energy which causes cosmic 
accelerat ion and constitutes 7 5% of the Universe' s energy 
density (|Albrecht et alj l200fj ; [Peacock et all 120061 ). There 
exists a wealth of potential models for dark energy. To dis- 
tinguish these models the determination of the dark en- 
ergy equation of state w has gained importance since some 
models can result in very different expansion histories. Cur- 



rent data can constrain the dark energy equation of state 
w to 10%, with the assumption of flatness, but a percent- 
age level sensitivity as well as redshift evolution information 
are required in order to understand the nature of dark en- 
ergy. Future cosmic shear surveys show exceptional poten- 
tial for constraining the dark energy equa tion of state w(z) 
(jAlbrecht et al. 2006; Pe acock et alj|2006h and have the ad- 
vantage of directly tra cing the dark matter distribution (see 
iHoekstra fc Jainll2008l for a review). 

In fact, cosmic shear surveys have the potential to 
constrain all sectors of our cosmological model. As shear 
measurements depend on the initial seeds of structure, it 
can be used to probe the slope and runn ing of the initial 
power spectrum (see e.g. iLiu et al]|2009l ) which is central 
to our understandin g of the inflationary model (see e.g. 
iHamann et al.l 120071 ) . Shear measurements have also been 
used to complement neutri n o constraints from particle 
phys ics (jTereno et alj 120091 ; llchiki. Takada. fe Takahashi 



l2009h and galaxy surveys (Tak ada. Komatsu. fc Futamasd 
2006), and future weak lensing surveys will provide bounds 
on the sum of the neutrino masse s, the number of massive 
neutrinos and the hierachy f kannesta d. Tu. fc Wonsll^OOd; 
Kitching. Taylor, fc Heavenj l2008bl ; |Pe Bernardis et al.l 
20091) . 
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The parameters that describe each sector are degener- 
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ate in the cosmic shear power spectrum, which means fixing 
parameters in one sector may results in anomalies being de- 
tected in another. We think all sectors should therefore be 
considered simultaneously. As each sector provides informa- 
tion about a different area of physics, we also argue they 
should also be considered on equal footing, i.e. have roughly 
the same number of parameters describing them, so that one 
sector is not favoured. Indeed, evidence for a departure from 
ACDM may come from any sector. 

In section 2, we describe the cosmological parameter set 
we consider, which includes the three sectors of dark energy, 
initial conditions and neutrinos. We introduce a new param- 
eter to include the second order running of the initial power 
spectrum. We also present the weak lensing tomography and 
the Fisher matrix forecast methods, and briefly discuss sys- 
tematic effects. In section 3, we present the weak lensing 
constraints we expect from future space based surveys with 
and without Planck priors and investigate the stability of 
the results as new parameters are added to the analysis. 
In section 4, we optimise such a survey using the Figure of 
Merit as well as constraints on all sectors. In section 5, we 
present our conclusions. 



2 METHOD 
2.1 Cosmology 

We start by describing the 12-parameter cosmological model 
which includes dark energy, dark matter (hot and cold) and 
initial conditions sectors. Throughout this paper, we work 
within a Friedmann-Robertson- Walker cosmology. Our cos- 
mological model contains baryonic matter, cold dark matter 
(CDM) and dark energy, to which we add massive neutrinos 
(i.e. hot dark matter - HDM). We also consider different 
parametrizations of the primordial power spectrum (defined 
in section |2~2"1 

We allow for a non-flat geometry by including a dark en- 
ergy density parameter Ode together with the total matter 
density O m , such that in general O m + Ode 1. The dy- 
namical dark energy equation of state parameter, w — p/p, 
is expressed as function of redshift and is parametrized 
by a first-order Taylor expansion in t he scale factor a 
jChevallier fc Polarski|[200ll ; lLindeijl2003l ): 

w(a) =w n + {a„ - a)w a , (1) 

where a — (1 + z) -1 . The pivot redshift corresponding to a n 
is the point at which w a and w n are uncorrelated. 
Our most general parameter space consists of: 

(i) Total matter density - O m (which includes baryonic 
matter, HDM and CDM) 

(ii) Baryonic matter density - 0;, 

(hi) Neutrinos (HDM) - m v (total mass), N v (number of 
massive neutrino species) 

(iv) Dark energy parameters - Ode, wq, w a 

(v) Hubble parameter - h 

(vi) Primordial power spectrum parameters - <rg (ampli- 
tude), n s (scalar spectral index), a (running scalar spectral 
index), f3 (defined in section |2~2")) 



We shall refer to 
'QCDM + a + p\ We 



ues 



based on the five-year WMAP result 



this fiducial cosmology as 
choose fiducial parameter val- 
ir WMAP results l|Dunklev et all 



I -MM H i similar to those used in iKitching etahl l|2008bh . The 
values are given in Table [1] 



2.2 Matter power spectrum 

The matter power spectrum is defined as: 
{5(k)S*(k')> = (2^) 3 «5?,(k-k')P(ft) 



and can be modelled by: 
P(k, z ) = ^A s k"° w+3 T 2 (k 



z) {d(o)J 



(2) 



(3) 



where A 3 is the normalisation parameter, T(k, z) is the 
transfer function and D(z) is the growth function. The pri- 
mordial spectral index is denoted by n s (k), and can depend 
on the scale k. 

In our cosmological model, the shape of the primor- 
dial power spectrum is of particular interest, since it may 
mimic some of the small-scale power damping effect of mas- 
sive neutrinos. In the concordance model, the primordial 
power s pectrum is generally parametrized by a pow er-law 
(see e.g. iKosowskv fc Turnerl fl995; Brid le et al.ll2003T) 



V x {k) = A s 



k 

k s 



(4) 



We parametrize the running of the spectral index by us- 
ing a second-order Taylor expansion of V x in log-log space, 
defining the running as a = dn s /d In k\k , so that the 
primordial power spectrum is now sc ale-dependent , with 
the scalar spectral in dex defined by (S pergel et alj 120031 ; 
lHannestad et alj|2002l ) 



n 3 (k)=n 3 (k ) + l^ 



In I — 

k 



(5) 



where fco is the pivot scale. We use a fiducial value of 
fco = 0.05Mpc -1 for the primordial power spectrum pivot 
scale. 

Although it is motivated by simplicity and standard 
slow-roll inflation theory, the second-order truncated Tay- 
lor expansion is limit ed and may lead to incorrect param- 
eter estimation (see lAbazaiian. Kadota. fc Stewar t 2005; 
iLeach fc Liddk]|2003h . In order to test this, we allow an ex- 
tra degree of freedom in the primordial power spectrum by 
adding a third-order term in the Taylor expansion, which we 
call p : 



1 ( k 

n s (k) = n s (k ) + -rain — 
2! \ko 

where f3 = d 2 n s / dln k 2 



1 , ( k 
+ 3! /31n U 



(6) 



We use the lEisenstein fc HvJ l|l999l ) analytical fitting 
formula for the time-dependent transfer function to calcu- 
late the linear power spectrum, which includes the contri- 
bution of baryonic matter, cold dark matter, dark energy 
and massive neutrinos, with the modification in the t ransfe r 
function su ggested by | Kiakotou. Elgar0v. fc Lahavl (2008). 
We use the [Smith et all (|2003l )" correction to calculate the 
non-linear power spectrum. The matter power spectrum is 
normalised using as, the root mean square amplitude of the 
density contrast inside an 8/i -1 Mpc sphere. 
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Table 1. Cosmological parameter sets used in our calculations. For each parameter set, the ticks (/) and crosses (x) indicate whether 
a parameter is allowed to vary or not, respectively. 



Parameters 


Wo 


W a 


^DE 


S2m 


Q b 


h 


0-8 




a 





m v 


Nu 


Fiducial values 


-0.95 





0.7 


0.3 


0.045 


0.7 


0.8 


1 








0.66 


3 


QCDM 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


X 


X 


X 


X 


QCDM + a 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


X 


X 


X 


QCDM + Q + /3 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


X 


X 


yQCDM 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


X 


X 


/ 


/ 


i/QCDM + a 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


X 


/ 


/ 


i/QCDM + a + P 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 



Following lEisenstein fc Hul (|l999h . we assume N„, the 
number of massive (non-relativistic) neutrino species, to be 
a continuous variable, as opposed to an integer. 

Neutrino oscillation experiments do not, at present, de- 
termine absolute neutrino mass scales, since they only mea- 
sure the difference in the squares of t he masses between 
neutrino mass eigenstates ijQuigdfeoOct ). Cosmological ob- 
servations, on the other hand, can constrain the neutrino 
mass fraction, and can distinguish bet ween different mass 
hierarchies f se e lElgarOv fc Lahavl [2005 for a review of the 
methods). The lEisenstein fc Hul transfer function assumes a 
total of three neutrino species (i.e. iV mass i css +N U = 3), with 
degenerate masses for the most massive eigenstates, i.e. if 
m„ is the total neutrino mass, then 



E 



mi = N v m,i 



(l) 



where rm is the same for all eigenstates. Thus, N v — 2 for 
the normal mass hierarchy, and N v — 1 for the inverted 
mass hierarchy, while N v — 3 corresponds to the case where 
all th ree neutrino species have the same mass (see Quige 
2008). The temperature of the relativistic neutrinos is as- 
sumed to be equal to (4/ll) 1//3 of the photon temperature 
l|Kolb fc Turnerlll990l ). 

Dark energy affects the matter power spectrum in three 
ways. Its density SIde changes the normalisation and k eq , 
the point at which the power spectrum turns over. Que and 
the dark energy equation of state parameter w change the 
growth factor at late times by changing the Hubble rate. In 
addition to this, for departures from a cosmological constant 
the shape of the matter power spectrum on large scales is 
affected through dark energy perturbations. In all our cal- 
culations, we only consider small deviations from wo — — 1, 
and so neglect dark energy perturbations, only considering 
the first two mechanisms. 

In comparing parameter constraints in different param- 
eter spaces, we shall use six parameter sets. We start with 
the simplest set (QCDM) to which we add neutrino and ad- 
ditional primordial power spectrum parameters. In all cases 
the central fiducial model (given in Table [1} is the same, but 
the number of parameters marginalized over varies. 



2.3 Weak lensing tomography 

The cosmological probes considered in this paper are tomo- 
graphic cosmic shear and the CMB. In weak lensing sur- 
veys, the observable is the convergence power spectrum. 
In our analysis, we calculate this quantity from the mat- 
ter power spectrum via the lensing efficiency function. Our 



Table 2. Fiducial parameters for the all-sky weak lensing survey 
considered. 



A s /sq degree 


20 000 


^■median 


0.9 


n s /arcmin 2 


35 


a z {z)/(l + z) 


0.025 


c e 


0.25 



convergence power spectrum therefore depends on the sur- 
vey geometry and on the matter power spectru m. We use 
the p ower spectrum tomography formalism by iHu fc Jam] 
(2004), with the background lensed galaxies divided into 10 
redshift bins. Cosmological models are then constrained by 
the power spectrum corresponding to the cross-correlations 
of shears within and between bins. The 3D power spectrum 
is proje c ted on to a 2D lensing correlation function using the 
iLimberl (|l953l ) equation: 



°t = I dz^ r W i (z)Wj{z)P(k = i/D A ,z), 



(8) 



where i, j denote different redshift bins. The weighting func- 
tion Wi(z) is defined by the lensing efficiency: 



Wi{z) 



3 Ho HqDol 



dz 



L>OS 



(9) 



where the angular diameter distance to the lens is Dol, the 
distance to the source is Dos, and the distance between the 
source and the lens is D LS (see lHu fc Jalnll2004l for details). 
Our multipole range is 10 < I < 5000. 

The galaxies are assumed to be distributed accord- 
ing to the following proba bility distribution function 
<|Smail. Ellis, fc Fitchettlll994h : 



P(z) — z a exp 



(10) 



where a = 2 and b = 1.5, and zo is determined by the media n 
redshift of the (see e.g lAmara fc Refregierll200?t ). 

Our survey geometry follows the parameters for a 'wide' 
all-sky survey with A s = 20000 sq degrees. The survey pa- 
rameters are shown in Table [2] The median redshift of 
the density distribution of galaxies is z me dian an d the ob- 
served number density of galaxies is n g . We include photo- 
metric redshift errors a z (z) and intrinsic noise in the ob- 
served ellipticity of galaxies <r e . We follow the definition 
cr~ = of , where a-y is the variance in the shear per galaxy 
fsee lBartelmann fc Schn eider 200l|). 
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2.4 Error forecast 

The predictions for cosmological parameter errors pre- 
sented in this paper use the Fisher matrix formalism. The 
Fisher matrix gives us the lower bound on the accuracy 
with which we can estimate model parameters from a 
given data set (Fisher 1935j; iTegmark. Taylor, fc Heavens! 
Il997l ; iKitching fe Amarall2009l ). In calculating forecast sur- 
vey errors, we are implicitly making assumptions about 
th e parameter set (see the d i scussi on of nested models 
m Heavens, Kitching. fc Verde! 120071 ). We want to know 
whether our constraints are robust against variations in the 
parameterisation of the cosmological model. This is of par- 
ticular importance when dark energy constraints are consid- 
ered, because of the degeneracies with other parameters (see 

Fig. mi. 

The forecast parameter precision is improved by com- 
bining independent experiments. Using this technique, joint 
constraints or error forecasts can be obtained by combining 
weak lensing with other observational techniques. 

The Fisher matr ix for the shear power spectrum is given 
by (|Hu fc Jainll2"00i ): 



F, 



/sky 22 



(2£ + l)A£ 



Tr [d^C^D^C^ 



(11) 



where the sum is over bands of multipole I of width A£, Tr 
is the trace, and / s k y is the fraction of sky covered by the 
survey. Equation [11] assumes the likelihod obeys a Gaussian 
distribution with zero mean. The observed power spectra for 
each pair i,j of redshift bins are written as the sum of the 
lensing and noise spectra: 



The derivative matrices are given by 

[D a r 



(12) 



(13) 



where p a is the vector of parameters in the theoretical 
model. 

In order to quantify the potential for a survey to con- 
strain dark energy parameters, we use the Figure of Merit, 
as de fined by the Dark Energy Task Force |Albrecht et al.l 
120061 ): 



FoM 



1 



Aw n Aw a 



(14) 



2.5 Planck priors 

In this article, together with lensing-only constraints, we 
also include joint lensing and CMB constraints. To com- 
bine constraints from different probes we add the respec- 
tive Fisher matrices: Fj i Ilt = Fi cns i ng + Fqmb- For our 
CMB priors, we use the forthcoming Planck mission as 
our su rvey. The Pl a nck F isher matrix is calculated fol- 
lowing iRassat et alj (|200ST ) . which estimates errors using 
information from the temperature and E mode polarisa- 
tion (i.e. TT, EE, and TE), with the help of the publicl y 
available camb code l|Lewis. Challinor. fc Lasenbvl |2000[) . 
We conservatively do not use information from B modes, 
and only use the 143 GHz channel, assuming other fre- 
quencies will be used for the foreground removal. This is 



conservative compared to other Planck priors i n the liter- 
ature. More details are given in IRassat et af] (|2008l . Ap- 
pendix B). The full parameter set for the Planck calculation 
is: {Ode, wo, w a , Q m , Q b , m„, iV„, h, a 8 , n s , a, t}. We use 
the same central values as for our weak lensing calculations, 
as described above, with a fiducial value for the reionisation 
optical depth r = 0.09 and subsequently marginalize over r. 
We consider neutrino parameters as Q v h 2 and N v and use 
a Jacobian to translate this into constraints on m v and N v . 



2.6 Systematic effects 

Weak lensing measurements are affected by systematic ef- 
fects which reduce the p recision on cosm olo gical parameters 
and introduce bias (see Rcfrcgicr 2003 and Schneider 2006 
for a review). In this section we discuss some of these effects. 

Intrinsic correlations can contaminate the lensing sig- 
nal. Solutions include using tomography or 3D lens- 
ing, which decouple the long-distance line-of-sight ef- 
fects from the the physical prox i mity of the galax- 
ies (see e.g. iKing fc Schneider! l2003l : IBridle fc King! 120071 : 
IKitching et alj |2008b). Using this approach, a nulling tech- 
nique for shear-intrin s ic ell i pticity has been proposed by 
iJoachimi fc Schneider! (|2008i . l2009h . 

Measurement systematics are due to PSF effec ts (see 
iKaiser et all Il995l ). The results presented in IBridle et alj 
|2009l) indicate that the required accuracy will be reached 
for the next generation of all-sky weak lensing surveys. 

Redshift distribution systematics, which lead to an 
uncertainty in the median galaxy redshift, cause an un- 
certainty in the amp litude of the matter power spectrum 
l|Hu fc Teg mark 1999). The problem of photometric redshift 
systematics can be met given a number of galaxies in the 
spectroscopic calibration sample of 10 4 — 10 5 (|Ma et al.ll2006l : 
lAmara fc Refregied[2007l '). 

Finally there are theoretical uncertainties on the shape 
of the matter power spectrum. The exist ing non-linear cor- 
rections to the matter power spectrum ([Peacock fc Dodda 
1 19961 : [Mai Il998t ) are only accurate to about 10% and dis- 
agre e with one ano ther to this level in the non-linear regime 
(see Hf utcrcr 2 00 ll ). Newer prescriptions such as those by 
ISmith et all (2003) offer more accurate predictions partic- 
ularly for non-ACDM cosmological models. The error in 
the non-linear part may still be si gnificant if effect of mas- 
sive n eutrinos is included (see e.g. ISaito. Takada. fc Taruval 
2008). Current semi-analytical models need to be im- 
proved to match the degree of statistical accuracy ex- 
pected for future weak lensing surveys. The solution is 
to run a suite of N -body ray-tracing simula t ions (see e.g. 
I White fc Valell2004l: iHuterer fc Takadal 120051 : iHilbert et all 
l2009l : ISato et al.ll2009l : lTevssier et al.ll2009f ). 



3 RESULTS 

In our Fisher matrix formalism, the error forecast on each 
parameter depends on the sensitivity of the weak lensing 
observation to changes in the matter power spectrum. In 
order to probe the effect of the different parameters on the 
matter power spectrum, we consider the fractional change in 
the non-linear matter power spectrum P(k), defined as the 
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change in P(k) with respect to the fiducial P(fe)ad, when 
one parameter at a time is varied from its fiducial value: 

Fractional change = P ^~ P ^ , (15) 
P(k)fid 

where A = 10% for all parameters in the ^QCDM + a + (3 
parameter set The power spectrum is normalised using as- 
The fractional change in P(k) at redshift z — is shown 
in Fig. [T] There are several features of interest in this plot, 
including the degeneracy between the parameters a, ft, m v 
and fi;, at small scales, as well as the degeneracy between 
wo, w a , f^DE and N v at large scales. The plot shows that 
the non-linear matter power spectra for the fiducial model 
and for the model with non-zero w a are almost completely 
degenerate at z = 0. This degeneracy is lifted as the redshift 
increases. 

Table [3] shows the marginalized errors for each pa- 
rameter in our six cosmological parameter sets. Joint lens- 
ing-)- Planck marginalized errors are shown in Table [4] When 
expected errors for n unknown parameters are calculated us- 
ing a Fisher matrix, we are implicitly setting the errors on 
any additional parameters to zero. We should therefore ex- 
pect QCDM to give us the best parameter constraints. In 
order to examine the variation in the marginalized errors 
with respect to the 8-parameter QCDM set, we define the 
fractional change in the marginalized error A for each pa- 
rameter as 

Fractional change in error = ^^ cx * — ^PQCdm^ 

Apqcdm 

where the subscripts cxt and qcdm denote the 'extended' 
hypothesis space and QCDM (our most restricted space) 
respectively. This quantity is shown in Fig. [2] for the eight 
parameters common to all the parameter sets, while Fig. [4] 
shows the fractional error change in the joint lensing+ Planck 
constraints. 



3.1 QCDM results 

The marginalized errors for the eight parameters in QCDM 
are shown in the second column of Table [3] With QCDM, 
all the sectors of our model are constrained well, even with 
our fiducial model containing massive neutrinos. Using our 
fiducial weak lensing survey with a QCDM parameter set, we 
obtain ~ 5% expected precision on wq. The joint errors on 
the dark energy parameters wo and w a are shown in Fig. [5] 
The FoM in this case is 130.99. With the addition of Planck 
priors, we find a significant improvement in the error bounds 
for the fioE, fim, fit, h, as, and n s . The improvement in the 
error bounds on wo and w a is smaller, with the FoM being 
increased by a factor of 2.75 (Tabled second column). 

3.2 Neutrino parameters 

In our calculations, we constrain the total neutrino mass and 
the number of massive species by measuring their effect on 
the lensing power spectrum via the matter power spectrum. 
This is sensitive to the neutrino fraction, r elated to the total 
neutrino mass by (|Elgar0v fc Lahavll2005l ) 

f - n " - 1 ( m »\ fly) 
3v fi m 94fi m h 2 VeV/ ' y ' 



The matter power spectrum parameterisation is also sensi- 
tive to_thejiuinber ofmassive neutrino species N u . 

iTereno et aL | <|2009h find a 3.3 eV upper bound for 
the total ne u trino mass, using CFHTLS-T0003 data, while 
llchiki et all (|20Qgh find an upper bound of 8.1 eV. Using 
our fiducial i/QCDM + a cosmology with neutrino parame- 
ter values of m v = 0.66 eV and N v — 3, our marginalized 
error forecast for m v is 1.20 eV, which gives us a la upper 
bound of 1.86 eV for the total neutrino mass (see Table [3J. 
With our joint lensing and Planck constraints (Table [1J, we 
obtain an error of 0.14 eV. 

The error on the dark energy density fioE is stable 
against the addition of massive neutrinos to the parame- 
ter set. For the equation of state parameters, we observe 
a degradation in the marginalized constraints. The FoM is 
consequently also degraded, as can be seen in Table [3] The 
top panel of Fig. [2] shows that the parameter most sensitive 
to the addition of neutrinos is Hubble parameter h, and to 
a lesser extent, wo and w a . In the latter case, this is due to 
a degeneracy with neutrinos in the observed effect on the 
growth function. 

Fig. [5]shows the joint la constraints on the dark energy 
parameters Wo and w a . The addition of massive neutrinos to 
the QCDM set produces a degradation on these constraints 
but does not significantly change the orientation of the el- 
lipse. This means that the pivot point a n remains almost 
unchanged. 

3.3 Constraints on the primordial power spectrum 

In this section, we find constraints for primordial power spec- 
trum parameters, using different parameterisations of the 
primordial spectral index. We also examine the variations in 
our error forecasts for other parameters when we vary our 
primordial power spectrum parameterisation. Weak lensing 
forecast co nstraints on t he ru nning spectral index have been 
studied bv llshak et ail (|2004 ). who also use combined lens- 
ing + current CMB constraints. They find that CMB con- 
straints are improved when weak lensing is added, especially 
for the parameters as, fi m , h and SIde- Their cosmological 
model does not include massive neutrinos, however. 

The parameterisation of the primordial power spectrum 
used here assumes a pivot scale feo at which the amplitude 
is defined. We find that the best constraints on n s in the 
i^QCDM + a set using our fiducial weak lensing survey are 
achieved with a pivot scale fco ~ IMpc -1 , which is larger 
than the value of 0.05 Mpc -1 , adopted in the rest of the 
paper. The results are shown in Fig. [3] This optimum pivot 
scale is shifted to ~ 0.3 Mpc -1 when the parameter /3 is 
added, and to ~ 0.1 Mpc -1 when a is added. 

In Table UJ we note that the addition of the parameter 
a has a small effect on the FoM, while adding a further pa- 
rameter P produces a larger degradation. The degradation 
in the fij, , fi m and as constraints is negligible against the ad- 
dition of a, while the parameters h and n s are most affected, 
as can be seen in Fig. [2] bottom panel. The addition of f3 
degrades the constraints on all these parameters, especially 
n 3 . 

With weak lensing only, we obtain tighter constraints 
on a than on n a with the QCDM + a and QCDM + a + /3 
parameter sets. This error hierarchy is reversed when Planck 
priors are added (Tables |3j and |4j fourth and fifth columns). 
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Figure 1. The fractional change in the non-linear matter power spectrum P(k,z = 0), obtained by varying each parameter in the 
i/QCDM + a + set by +10% from its fiducial value. 



Table 3. Predicted marginalized parameter errors for weak lensing alone. We show results using different cosmological parameter sets. 
The second column shows the results for our most restricted parameter set QCDM. In the third column, we add massive neutrinos. 
Primordial power spectrum parameters are added in the fourth and fifth columns. In the sixth we add neutrinos and a running of the 
primordial spectral index. The seventh column shows our most extended set ^QCDM + a + /?. We also show the DETF Figure of Merit 
for each set. 



Parameter 


QCDM 


i/QCDM 


QCDM 
+a 


QCDM 

+a + /3 


v QCDM 
+a 


z^QCDM 
+a + P 


wo 
w a 
^DE 


0.05633 
0.19297 
0.05214 


0.06443 
0.23674 
0.05841 


0.05740 
0.21567 
0.05287 


0.06583 
0.24988 
0.05297 


0.08099 
0.32904 
0.05842 


0.09608 
0.48144 
0.05856 


n b 

m„/eV 


0.00731 
0.02411 


0.00742 
0.02558 
1.10229 
3.27380 


0.00731 
0.02544 


0.00752 
0.02981 


0.00749 
0.03200 
1.19614 
3.81643 


0.00756 
0.03201 

I. 51694 

II. 12214 


h 


0.11337 


0.23176 


0.18660 


0.24691 


0.41999 


0.53253 


n s 
a 
P 


0.01184 
0.02904 


0.01230 
0.03038 


0.01185 
0.08662 
0.04378 


0.01268 
0.11158 
0.05661 
0.02574 


0.01307 
0.11969 
0.06556 


0.01319 
0.12003 
0.07137 
0.08479 


FoM 


130.99 


79.69 


114.86 


97.59 


56.36 


38.02 



Examining Fig. [S] we observe that the addition of pri- 
mordial power spectrum parameters produces a small degra- 
dation in the joint (wo, w a ) errors and has little effect on the 
orientation of the ellipses. 



3.4 Combined neutrino and primordial power 
spectrum parameters 

We also investigate the effect adding both neutrinos and 
primordial power spectrum parameters (the sets i/QCDM+a 
and fQCDM + a + /5). We note that the effect on the FoM is 
more significant than with neutrinos or a and /3 alone (Table 
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Figure 2. The fractional change in the marginalized error for each parameter with respect to the QCDM set, using a lensing Fisher 
matrix calculation. In the top panel, we add neutrino parameters before adding degrees of freedom in the primordial power spectrum. 
In the bottom panel we add primordial power spectrum parameters before adding neutrinos. In each case, the most general parameter 
space, shown by the red bars, is QCDM + u + a + 0. The la marginalized errors are calculated using Fisher analysis for our all-sky 
fiducial weak lensing survey. 
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Figure 3. Marginalized error A on primordial power spectrum 
parameters against pivot scale feo- The calculation was carried 
out for two parameter sets: f QCDM + a (shown in black) and 
i^QCDM + a + f3 (shown in red) using our all-sky fiducial weak 
lensing survey. We show the marginalized errors for the parame- 
ters n s (solid line), a (dashed line) and /3 (dotted line). 

[3]). With the full extended model, the effect is especially 
noticeable on N v and /3, showing that there are significant 
degeneracies between the effect of neutrinos on the matter 
power spectrum and the effect a scale-dependent primordial 
power spectrum with several degrees of freedom. 

Fig. [2] (red bars) shows that the greatest degradation 
in constraints with respect to QCDM occurs in the param- 
eters w a , h, and n B . There is an additional degeneracy in 
the matter power spectrum between the small-scale power- 
suppression effect of massive neutrinos and the form of the 
primordial power spectrum for certain values of the primor- 
dial spectral index. 

With weak lensing only, we obtain tighter constraints on 
a than on n 3 , which is in a gree ment with the resul ts obtained 
bv lKitching et all (|2008al ) and lfshak et~ai1 (|2004l ). Although 
the precision for both parameters is degraded when neutri- 
nos are added, this error hierarchy is preserved, even when 
CMB constraints are added (see section 1531 below) . 

3.5 Joint lensing and CMB results 

The addition of Planck priors has a significant effect on pa- 
rameter constraints. The FoM is improved by a factor of 6 for 
the z/QCDM + a + ft model (Table @}, and we obtain better 
constraints for all parameters, especially Qb, ^2db (related 
to the geometry of the universe), h and n 3 . Adding CMB 
priors also lifts the degeneracy between some parameters, 
so the extension of the parameter set does not significantly 
degrade the error bars. This can be seen in Fig. fj] where 
Sib, h and n s , which are well-constrained by Planck, are now 
hardly affected by the addition of extra parameters, ft can 
be seen from Table fj] that we obtain better constraints on a 
than on n s with the addition of CMB priors, reversing the 
error hierarchy obtained with lensing only. Moreover, the 
addition of neutrino parameters does not significantly affect 
the precision on n s . 

With combined lensing+P/ancfc calculations, we obtain 




-1.10 -1.05 -1.00 -0.95 -0.90 -0.85 -0.80 



Wo 

Figure 5. Joint la constraints in the (wo,w a ) plane from our 
Fisher matrix calculation with four different parameter sets. 
QCDM, i/QCDM, QCDM + a-l-^ and uQCDM + a + /3 are shown 
in grey, red, blue and green respectively. The solid ellipses show 
the constraints using lensing only, while the dashed ellipses show 
the constraints with the addition of Planck priors from our CMB 
Fisher matrix calculation. 

an improvement in the joint (wq, w a ) constraints (Fig. [5}. 
The constraints are robust against the addition of neutrino 
parameters and the primordial power spectrum parameters 
a and (3. 



4 DEPENDENCE OF PRECISION ON 
SURVEY DESIGN 

In this section we study the effect of the survey design 
on our error forecasts. The optimisation for an all-sky to- 
mographic weak lensing survey has been investigated by 
l|Amara fc Refregierll2007i ). who use the dark energy FoM 
as the optimisation benchmark. Our survey configuration is 
determined by the parameters: the area A s , the median red- 
shift z m , and the observed number density of galaxies n g . 
The lensing correlation function is additionally defined by 
the range of multipoles over which it is measured. In order 
to investigate the dependence of the marginalized errors in 
our i^QCDM + a + ft parameter set on the survey param- 
eters, we calculate the lensing Fisher matrix while varying 
one survey parameter at a time. In Fig. [6] we show the rel- 
ative marginalized error, defined as Ap/\p\ for different val- 
ues of the median redshift z m and the maximum multipole 
^max- We note that the scaling with z m is similar for all 
the parameters in our three sectors of interest, with N v and 
P showing a stronger dependence on z m . We find that all 
parameters have a roughly similar scaling with ^ ma x, and 
that the greatest gain in precision is observed in the range 
10 2 < ^ max < f0 4 . 

We also carried out the same calculation for various 
values of the survey area A s and the galaxy count n g , finding 
a linear scaling between the parameter precision and these 
two survey parameters. Our results show that the survey 
area has the great est effect on parame t er pre cision, and are 
in agreement with I Amara fc Refregierl (|2007l ) . 

Our results show that the optimum survey strategy 
holds not just for the joint dark energy parameters (wo,w a ), 
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Figure 4. Fractional change in marginalized errors with respect to QCDM, using a joint lensing+ Planck Fisher matrix. As in Fig. [2] in 
the top panel we add neutrino parameters before adding the parameter a, and in the bottom panel we consider different parameterisations 
of the primordial power spectrum without adding neutrinos to our parameter space. 
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Table 4. Predicted marginalized parameter errors for for cosmic shear combined with Planck priors. 



Parameter 


QCDM 


uQCBM 


QCDM 


QCDM 


v QCDM 


i/QCDM 








+Q 


+a + P 


+a 


+OL + P 


■wo 


0.04942 


0.04984 


0.04943 


0.05055 


0.04987 


0.05142 




0.17943 


0.18231 


0.17946 


18260 


18275 


0.18482 


^DE 


0.00644 


0.00661 


0.00721 


0.00722 


0.00730 


0.00730 




0.00389 


0.00391 


0.00391 


0.00391 


0.00393 


0.00393 


n b 


0.00091 


0.00119 


0.00101 


0.00101 


0.00128 


0.00128 


nrtv /eV 




0.14172 






0.14172 


0.14176 


N v 




0.11694 






0.11821 


0.11924 


h 


0.00599 


0.01360 


0.00625 


0.00625 


0.01381 


0.01382 




0.00461 


0.00491 


0.00467 


0.00470 


0.00492 


0.00501 


n 3 


0.00332 


0.00549 


0.00356 


0.00360 


0.00557 


0.00563 


a 






0.00515 


0.00519 


0.00545 


0.00545 


P 








0.01779 




0.01834 


FoM 


357.12 


258.40 


357.01 


348.70 


251.51 


240.59 



but also for the other parameters in the cosmological model. 
A different survey design would therefore lead to a reseat- 
ing of the marginalized errors shown in Table |3j but would 
not significantly modify the results shown in Fig. [2] This 
indicates that the effect of extending the hypothesis space 
is independent of the survey design. 



5 CONCLUSION 

The main aim of this paper was to put three different sectors 
of the cosmological model: dark energy, dark matter (cold 
and neutrinos) and initial conditions on equal footing while 
forecasting constraints for a future weak lensing survey. To 
do this we introduce a new parameter P, which models the 
second order running spectral index. We have forecast er- 
rors for an all-sky tomographic weak lensing survey, and for 
weak lensing+ Planck, using different cosmological parame- 
ter sets, studying the effect of the addition of parameters 
in the model. We have shown that error forecasts for some 
parameters are stable against changes in the parameter set 
(Table [3} , and that degeneracies between the dark energy 
parameters Wq and w a are not significantly affected by the 
addition of parameters (Fig. [SJ. We have investigated the 
shift in the optimal pivot scale for primordial power spec- 
trum constraints with the addition of extra parameters (Fig. 
0. 

We have shown that parameter constraints are never- 
theless dependent on the parameterisation of the primordial 
power spectrum. With the addition of CMB priors, we have 
shown that we can obtain improved constraints in this sec- 
tor, which reduces the dependency on constraints on the 
parameter set. 

We have obtained predicted constraints for both the 
total neutrino mass and the number of massive neutrino 
specie s. This article has a resonance with lDe Bernardis et al.l 
(2009), in which it is shown that the neutrino mass hierar- 
chy can be constrained using cosmic shear, using a more 
general parameterisation of the neutrino mass splitting. We 
find that for the parameters that are common between the 
two articles there is an agreement between the predicted 



errors, despite the slightly different parameter sets and as- 
sumptions. The results ob tained here also mirror those in 
IZunckel fc Ferreiral l|2007t) . in which it is found that the 
neutrino mass constraints are degraded when the hypoth- 
esis space is enlarged. This is due to degeneracies between 
the neutrino mass and other parameters. Adding priors only 
tightens the constraints when the additional information 
comes from independent experiments, which reduces the 
freedom in the degenerate parameters. We have used in- 
formation from the CMB, which constrains the primordial 
power spectrum particularly well, resulting in improved con- 
straints on our neutrino parameters. 

In the neutrino sector, parameter constraints would be 
improved by a hierarchical parameterisation in which differ- 
ent neutrino species have non-degenerate masses. This would 
model more accurately the process whereby each massive 
species becomes non-relativistic at a different redshift. While 
the different transition redshifts have only a very small effect 
on the CMB anisotropy power spectrum, the effect is non- 
negligible in future cosmic shear experiments which measure 
the matter power spectrum to a sufficient accuracy to dis- 
criminate between different mass hierarchies. 

This article concludes that a future all-sky weak lensing 
survey with CMB priors provides robust constraints on dark 
energy parameters and can simultaneously provide strong 
constraints on all parameters. The results presented here 
show that error forecasts from our weak lensing survey are 
stable against the addition of parameters to the fiducial 
model, and that this stability is improved by adding CMB 
priors. 
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Figure 6. Dependence of cosmological parameter precision on the 
survey design. In the top panel, we vary the median redshift z m , 
while keeping other survey parameters constant. In the bottom 
panel, we vary the maximum multipole £ max , with the minimum 
multipole kept constant at i = 10. Each plot shows the relative 
marginalized error, defined as Ap/|p| for the twelve parameters in 
i^QCDM + a + p. The thick black line shows the (wo, w a ) Figure 
of Merit. 
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